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ABSTRACT 

We present the results of global VLBI observations of the BL Lac object PKS 0003—066 
at 2.3 and 8.6 GHz, over a 10 year period using data from new and archived VLBA, geodetic 
VLBI, and global VLBI (including the use of new disk-based recording and software correlation 
systems) observations. Inter-continental baselines resolve this source into a compact core and 
three jet components that move away from the core with proper motions ranging between '--^O.l 
mas/yr for the outer jet component and ~1.4 mas/yr for the fastest inner jet component, i.e. 
apparent speeds of ^^2.0 c to ^^21. 6 c respectively. This result is in contrast to previous studies of 
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this source which did not reveal significant jet motions from Hmited data sets. We find that the 
fast inner jet components catch up to the more slowly moving outer jet components and interact 
strongly, causing local brightening of the jet in the interaction. We document the body of VLBI 
data for this source, as well as other supporting radio observations, as a basis for future work 
that will investigate the underlying physics of the pc-scale jet in PKS 0003—066. 

Subject headings: BL Lacertae object: individual (PKS 0003-066, NRAO 005, J0006-0623) - galaxies: 
jets - radio continuum: galaxies - techniques: interferometric 



1. Introduction 

BL Lac objects are extragalactic sources char- 
acterised by weak or undetectable optical line 
emission and strong and variable flux and polar- 
ization in a wide range of wavebands from op- 
tical through radio (for review, e.g. KoUgaart 
1994). In the radio band, very long baseline inter- 
feromctry (VLBI) images show that the pc-scalc 
structures of BL Lac objects are dominated by 
compact, flat-spectrum cores from which (appar- 
ently) one-sided jets emerge, often showing knots 
of enhanced brightness with apparent superlumi- 
nal speeds (e.g Gabuzda, Pushkarev & Cawthorne 
2000). These properties are explained by Doppler 
boosted synchrotron radiation from a relativistic 
jet closely aligned to our line of sight (Blandford & 
Recs 1974). In the framework of the unified model 
of radio-loud active galactic nuclei (AGNs) (see 
e.g., Urry & Padovani 1995), the observed differ- 
ences in the properties of quasars, radio galaxies, 
and BL Lac objects are due to the orientation of 
the relativistic beam and an obscuring torus with 
respect to our line of sight. 

The BL Lac object PKS 0003-066 (NRAO 005, 
J0006-0623) at a redshift of 2; = 0.347 (Stickel, 
Fried & Kiihr 1989) is a bright compact radio 
source and has been observed with VLBI for many 
projects (e.g. Fey & Chariot 1997). The core 
is unresolved at arcsecond resolution (Ulvcstad et 
al. 1981, Perley 1982). Gabuzda, Pushkarev & 
Cawthorne (2000) observed PKS 0003-066 with 
a global VLBI array and showed that the pc-scale 
structure is resolved, with the major components 
strongly polarised, but couldn't identify any su- 
pcrluminal motions from comparison to previous 
VLBI observations. Kellermann et al. (2004) 
analysed the results of a 15 GHz (2 cm) multi- 
epoch Very Long Baseline Array (VLBA) program 
begun in 1994, to study the jets in 110 quasars and 
active galaxies. PKS 0003—066 is one of 18 BL Lac 
objects in that sample, Kellermann et al. (2004) 



identifing three jet components with no significant 
proper motions. Total flux density monitoring of 
this source with the ATCA at 1.4, 2.5, 4.8 and 8.6 
GHz shows rapid variability over a 2 - 3 year pe- 
riod, with a higher degree of variability at higher 
frequencies (Tingay et al. 2003), suggesting vari- 
ability of the pc-scale structures. 

To investigate this possibility, we have made 
new global VLBI observations of PKS 0003-066 
at 2.3 GHz. We have also analysed 25 other VLBI 
observations at 2.3 and 8.6 GHz from a 10 year se- 
ries of observations to study the flux density and 
structural variability of the source in more detail 
than has previously been possible. In this pa- 
per we describe these results and a comparison to 
the results of previous multi-frequency and multi- 
epoch VLBI experiments. In contrast to the re- 
sults in Kellermann et al. (2004), we find evidence 
for fast jet components ejected from the core which 
interact strongly with slower jet components fur- 
ther from the core and cause local jet brightening. 

In section 2.1, we describe new VLBI observa- 
tions at 2.3 GHz using a global VLBI array, new 
disk-based recording techniques, and software cor- 
relation techniques. In section 2.2 we summarize 
the results obtained from archival VLBI data. In 
section 3 we discuss the pc-scale structure and evo- 
lution of this source. 

Throughout this paper we use the following cos- 
mological parameters as used by Kellermann et 
al. (2004): Hq = 70 kms-iMpc-\ = 0.3, 
flA = 0.7. The corresponding luminosity distance 
is 1.939 Gpc at the redshift of PKS 0003-066, 
z = 0.347, and the linear scale is therefore 4.91 
pc/mas. A proper motion of 1 mas/yr corresponds 
to an apparent speed of 16c. 
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2. Observations and Results 

2.1. Global VLBI observations at 2.3 GHz 

2.1.1. Global disk-based VLBI observations 

PKS 0003—066 was observed with a global ar- 
ray of antennas on 2004 April 18, using 5 Aus- 
tralian telescopes, one South African telescope 
(Hartebeeestoek Radio Astronomy Observatory), 
and one Japanese telescope (Kashinia), as sum- 
marised in Table 1, for a period of 20 minutes. 
This experiment was designed as a test observa- 
tion for the software correlator described in section 

2.1.2. The observational setup utilised a number of 
different disk-based recording systems (as briefly 
summarised in Tabic 1) to record the Nyquist- 
samplcd and digitised telescope output as a func- 
tion of time. At the antennas that utilised the 
LBA Data Recorder (LBADR) ^ recording sys- 
tem (West 2004), an aggregate data rate of 256 
Mbps was recorded, corresponding to four 16 MHz 
bands, each sampled at 32 MHz and 2-bit digi- 
tised. The four 16 MHz bands were arranged with 
two bands at right circular polarisation (RCP) 
(2274 - 2290 MHz) and two bands at left circu- 
lar polarisation (LCP) (2274 - 2290 MHz). Dupli- 
cate bands were recorded at RCP and LCP for 
redundancy. At the antennas that utilised the 
Mark5 recording system, an aggregate data rate 
of 128 Mbps was recorded, in matching 16 MHz 
bands. At the Kashima antenna, which utilized 
the Japanese K5 recording system (Koyama et 
el. 2003); an aggregate data rate of 1 Gbps was 
recorded. The 2274-2306 MHz frequency range of 
the RCP signal was converted to a 0-32MHz base- 
band signal, sampled at the rate of 512 MHz and 
then digitized with 2-bit precision. The 1 Gbps 
data was digitally processed to obtain the time se- 
ries data corresponding to a single 16 MHz band 
sampled at 32 MHz using 2-bit digitization. This 
digitally filtered 16 MHz band ranged from 2274 
- 2290 MHz, corresponding to one of the 16 MHz 
bands recorded using the LBADR and Mark5 sys- 
tems at the other antennas. An account of the 



^The LBADR recorders use the VLBI Standard Interface 
I/O Board (VSIB) and the VLBI Standard Interface Con- 
verter(VSIC) cards developed by the Metsahovi Radio Ob- 
servatory (Ritakari & Mujunen 2002), along with Apjjlc 
XRaid units for data storage, except Ceduna which used 
internally raided drives for this experiment. For further 
details, see Phillips et al. (2006, in preparation) 



software used to perform the digital filtering can 
be found in Takeuchi (2004). 

This experiment illustrates the power and flex- 
ibility of disk-based recording for VLBI. Tradi- 
tional tape-based VLBI requires considerable ef- 
fort to translate different recorder systems prior 
to correlation. Our software correlator is capable 
of reading different disk-based recorder formats 
without conversion, improving the compatibility 
of antennas with different recording systems. The 
data recorded at each antenna were shipped to the 
Swinburne University of Technology, for correla- 
tion using the software correlator described in sec- 
tion 2.1.2. Since the Kashima antenna only pro- 
duced a single 16 MHz band (RCP) using the soft- 
ware digital filter, only this band was correlated 
with a view to full data reduction and imaging. 

2.1.2. Description of correlation in software 

A software correlator for VLBI has been devel- 
oped using the Beowulf cluster at the Swinburne 
University of Technology Centre for Astrophysics 
and Supercomputing. The Beowulf cluster con- 
sists of a large number (~300) of Pentium based 
PCs (for detail, see http://astronomy.swin.edu.au/supercomputing/). 

The software correlator is an XF style correla- 
tor ^. Data recordc^d to hard disk at each individ- 
ual telescope are mounted on the supercomputer 
and the data for each antenna pair (baseline) is 
transfered to a processing node for correlation, in 
segments of up to 25 ms at a time. High level con- 
trol of the data flow is provided under the stan- 
dard Message Passing Interface (MPI). The data 
to be transfered is selected from the data on disk 
according to the geometric delay model and an- 
tenna clock models at the selected time in the ob- 
servation; the geometric delay model is based on 
an imaginary telescope located at the centre of the 
Earth, allowing the two data streams to be aligned 
to +/- 1 sample accuracy (the fractional sample 
error is retained for use in post-correlation correc- 
tions). The geometric delay model is generated 
using the CALC package (Himwich 1988, see also 
http: / / gemini.gsfc.nasa.gov/solve/solve.html) . 

Once the selected data from two antennas (a 



^This is a proto-type software correlator (West 2004). The 
correlator has more recently been superceded with a more 
capable FX software correlator (Deller et al. 2006 in prepa- 
ration) . 
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single baseline) are transferee! to a processing node 
for correlation, the data are unpacked into 32-bit 
floating point numbers. The fringe rotation func- 
tion derived from the predicted delay between the 
two antennas is then formed and applied to one 
of the two data streams. The two data streams 
are then cross correlated to form the lag spectrum 
(expressed in terms of the correlator coefficient 
at each lag), and the lag spectrum Fourier trans- 
formed to give the cross power spectrum. In order 
to recover the full signal to noise (e.g. D'Addario 
1989), a second fringe rotation function is gener- 
ated, 7r/2 out of phase from the first fringe rotation 
function. The correlation is repeated using the 
second fringe rotation function. The two resulting 
cross power spectra are then coherently averaged, 
after correction for the 7r/2 phase shift, to give 
the final cross power spectrum. A fractional sam- 
ple error correction is applied to the cross power 
spectrum, to correct for the phase slope across 
the band due to the limited accuracy in the ini- 
tial alignment of the two data streams. Finally, 
the measured system temperatures and gains from 
each antenna are applied to calibrate the visibility 
amplitudes into Janskys (see Table 1). 

The correlated data are reported back to a mas- 
ter node for collation and are exported to the FITS 
format, for further reduction in standard process- 
ing packages. This correlation scheme closely fol- 
lows that of the Australian Long Baseline Array 
(LBA) correlator (Wilson et al. 1992; Roberts 
et al. 1997). An account of a prototype ver- 
sion of the software correlator can be found in 
West (2004). This new system for global VLBI is 
available to all users through the Australia Tele- 
scope National Facility (ATNF) (for detail, see 
htt p : / / www . at nf . csiro.au/ vlbi / ) . 

2.1.3. Post- correlation data reduction and imag- 
ing analysis 

Following correlation, the data were imported 
into the AIPS package for fringe- fitting and phase 
calibration, then into the DIFMAP package (Shep- 
herd 1997) for editing, imaging, and evaluation 
of the source structure using modelfitting in the 
{u,v) plane. The {u,v) coverage for the obser- 
vation is thus shown in Fig. 1. The coverage is 
reasonably good for a 20 minute observation (es- 
sentially a single snapshot) due to the long east- 
west baselines from Australia to South Africa, 



the long north-south baselines from Australia to 
Japan, and the longest baseline from Japan to 
South Africa, which has large components in both 
the east-west and north-south directions. Fig. 2 
shows the visibility amplitude as a function of dis- 
tance in the (u, v) plane for the observation, pro- 
jected onto the u axis. Structure in the source is 
immediately apparent from the variations in the 
visibility amplitude. 

The data were model-fitted in the (m, v) plane 
using a three component model (Table 2) and self- 
calibrated in amplitude and phase. Amplitude 
corrections of less than 5% were derived from this 
procedure, showing that the a priori calibration 
(Table 1) was of a reasonable quality. The image 
resulting from the modelfit is shown in Fig. 3. 
As shown in section 2.2 and discussed in section 
3, the structure of PKS 0003-066 derived from 
these observations agrees well with the structure 
found from other VLBI array /correlator combi- 
nations such as VLBA-)-geodetic antennas at 2.3 
and 8.6 GHz (Fey & Chariot 1997, sec also be- 
low), the VLBA at 5 GHz (Fomalont et al. 2000, 
Gabuzda, Pushkarev & Cawthorne 2000), and 
the VLBA at 15 GHz (Kellermann et al. 1998). 
The core component is marginally resolved on 
the Australia-Japan-South Africa baselines, which 
are the longest baselines yet used to observe this 
source, apart from the VSOP observation for this 
source where three ground telescopes were used for 
two hours in conjunction with the HALCA space 
telescope (Scott et al. 2004) . The brightness tem- 
perature for the core component, Tb, is derived to 
be 5.3xlO^^K from Table 2, being consistent with 
the lower limits of > 4.8 x lO^^K derived by Scott 
et al. (2004) from the VSOP observations (see also 
Horiuchi et al. 2004) and > 4.60 x 10"K derived 
by Kovalev et al. (2005) from a VLBA 15 GHz 
observation. 

2.2. Analysis of VLB A/geodetic archive 
data 

Additional to the data described above, to 
study the structural variability of PKS 0003—066, 
we have used data from the Radio Reference 
Frame Image Database (RRFID) of the US Naval 
Observatory (USNO), at 2.3 and 8.6 GHz (Fey & 
Chariot 1997). Observations were made using an 
array consisting of the 10 antennas of the Very 
Long Baseline Array (VLBA) of the National Ra- 
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dio Astronomical Observatory (NRAO) along with 
geodetic antennas. In total 13 epochs of data are 
available from 1995 October to 2002 January (see 
Tabic 3). Each observation was undertaken with 
the dual-frequency VLBA receivers. The sources 
were observed using short duration (~ 3 minutes) 
"snapshots" over a number of different hour an- 
gles. We analysed and imaged using DIFMAP the 
original datasets for PKS 0003—066 kindly pro- 
vided to us by Alan Fey (private communication). 
As images from the RRFID are already published 
^, we do not reproduce them in full here. Full re- 
sults of our model-fitting analysis of the data are 
presented, however, in Table 4 for observations 
Usted in Table 3. 

We have also obtained two datasets from the 
VLBA online archive for PKS 0003-066. One is a 
part of a series of coordinated RDV (Research and 
Development VLBA) astrometric/geodetic experi- 
ments, numbering in total 9 epochs from 2002 July 
to 2004 July. These experiments use the full 10- 
station VLBA plus up to 10 Mark 4 geodetic sta- 
tions currently capable of recording VLBA modes. 
These experiments were coordinated by the geode- 
tic VLBI programs of three agencies: USNO, 
NASA, and NRAO to monitor source structure 
and to determine a high accuracy terrestrial refer- 
ence frame (Gordon 2000). 

Another dataset consists of three epoch obser- 
vations at 2.3 and 8.6 GHz with the VLBA only 
(obs. code BP 118) in May, June, and July 2004. 
These archival data are exported to AIPS and 
DIFMAP for imaging using standard procedures. 
Since these images from the VLBA archive data 
have not been previously presented, in Figs. 4, 5 
and 6 we show images obtained from our analysis 
of the 9 epoch RDV data and the 3 epoch BP 118 
VLBA data. Parameters of the images are given 
in Table 5. 

We thus have a total of 25 epochs of data at 2.3 
and 8.6 GHz in addition to the global VLBI obser- 
vation at 2.3 GHz described in Section 2.1. Fig- 
ures 7 and 8 show subsets of the RRFID and RDV 
images as montages at both 2.3 and 8.6 GHz over 
a 10 year period. The image from year 2003.35 is 
from the RDV data, images at all other epochs are 



^AU images are available from the USNO-RRFID web- 
site (see http://rorf.usno.navy.mil/RRFID/) and the 
NASA/IPAC Extragalactic Database (NED) 



from the RRFID data. 

From a model- fit analysis of the 8.6 GHz data 
we identified 3 components in the pc-scale jet (as 
did Kellermann et al. 2004) that are persistently 
recognisable until 1998, in addition to the core. 
Then a new, 4th, component emerged from the 
core in late 2002. Table 4 lists the source models 
we obtained at each epoch. Fig. 9 shows the an- 
gular distance of each component from the core as 
a function of time. Fig. 10 shows the offsets of the 
jet components in right ascension and declination 
from the core position, as a function of time. Fi- 
nally, the flux density variabihty of PKS 0003-066 
is shown in Fig. 11 for 8.6 GHz and Fig. 12 for 
2.3 GHz, using the all VLBI data analysed, decom- 
posing the source into its different pc-scale com- 
ponents. 

An interpretation of the flux density monitoring 
observations such as the Australia Telescope Com- 
pact Array (ATCA) (Tingay et al. 2003) and the 
University of Michigan 26m (UMRAO) database 
(see e.g. Aller et al. 1985 for description of the 
monitoring program) will appear elsewhere (Hori- 
uchi & Tingay 2006, in preparation). 

3. Discussion 

3.1. Proper motions and variability of the 

jet components 

From their analysis of 5 epochs of data over 
5.5 years (from July 1995 to January 2001) at 15 
GHz, Kellermann et al. (2004) fotmd no siginifi- 
cant proper motions for any pc-scale jet compo- 
nents in PKS 0003—066. In contrast, from an 
examination of the montages in Figs. 6 and 7, 
substantial evolution of the source structure is ap- 
parent when considering a large amount of data 
over a ten year period. We fit the radial dis- 
tances from the core versus time for each com- 
ponent of the 8 GHz images with the linear mo- 
tion that minimizes the (Fig- 9 (a)). The 
uncertainty of the proper motion, defined by a 
significance of > 3a, is estimated using the 
statistics that includes both scattering of the data 
points around the proper motion-model and un- 
certainties of individual observations (one quarter 
of the beam sizes are adopted). We show that 
the inner jet components, denoted C2 and C4 (see 
Fig. 9 (a)), have proper motions of 0.58±0.02 
and 1.35±0.05 mas/yr, corresponding to apparent 
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speeds of 9.3±0.3 c and 21.6±0.8 c. These inner 
components are much faster than the outer com- 
ponent, CI, which is moving at 0.12±0.01 mas/yr 
or approximately 2.0±0.2 c (Table 8). 

In Fig. 9 (b), the 15 GHz model fit components 
by Kellermann et al. (2004) based on their image 
analysis (no error estimates on their component 
position were provided) are superposed on the 8.6 
GHz fits for the linear motions together with the 
2.3 GHz model components. The 2.3 GHz compo- 
nents follow the similar trend of proper motions as 
seen at 8.6 GHz. Although there is a hint of shifts 
within one milliarcsecond or so between the paths 
at 2.3 and 8.6 GHz, the positional acuracy of 2.3 
GHz data, typically about 0.5 mas, is comparable 
to this. 

The reason why Kellermann et al. (2004) didn't 
identify such high proper motions may be at- 
tributed to the misidentification of the compo- 
nents due to the large gaps between observations. 
We have data points almost every two months 
from 1997 to 1998, hence we clearly see the fast 
proper motion of C2 contimiOTisly over that pe- 
riod. The fact that the results of proper motion 
studies of core-dominated radio sources can be so 
affected by the sampling frequency and overall pe- 
riod of monitoring is well known (e.g. Gomez et al. 
2001). Moreover, Fig. 9 (a),(b) indicates that our 
data point for C2 at the epoch of 16 January 2002 
(MJD 52291) is consistent with those for the com- 
ponent C of Kellermann et al. (2004) around that 
epoch. Also our images are at a lower frequency, 
hence steep spectrum components such as CI, C2 
and C4 are brighter (see Fig 12(b)) and easier to 
identify in spite of the lower spatial resolution. In 
contrast, at 15 GHz some jet components are re- 
solved and are weaker, making a consistent identi- 
fication of the overall jet structure more difficult. 
Note also that Kellermann et al. (2004) derived 
positions in the image plane while we derived po- 
sitions in the UV plane. This may lead to some 
differences in the derived positions for the compo- 
nents. 

3.1.1. Kinematics of the fast jet components 

Adopting the relation between the intrinsic jet 
speed (3 = v/c < 1, the apparent speed /3app) and 



the viewing angle of the jet to our line of sight 6 

^"PP ~ l-/Jcos^' 

we find that the jet inclination that maximises the 
apparent speed is ^max = 5.3°, using the fastest 
apparent motion, /3app = 21.6, of C4 or 12.3° for 
the 9.3c apparent motion of C2. One can then 
derive a minimum Lorentz factor, 7niin = 21.6 for 
C4, or 9.35 for C2, which corresponds to a Doppler 

factor of Smin = 7min(l " /3min COS 6>max )~^ =8. 63 

for C4 and 3.74 for C2, where /J^in = \/l~7^ = 
0.9989 for C4 or 0.9942 for C2. For smaller view- 
ing angles the Doppler factor approaches (5 ~ 27 
as 6* ~ 0, which yields (fmax = 43.2 for C4 and 18.7 
for C2. Although Kellermann et al. (2004) didn't 
measure superluminal motions in PKS 0003—066, 
in thier sample of 110 sources for the VLB A 2cm 
suevey they measued superluminal motions rang- 
ing up to l3app ~ 25. Our results suggest that the 
superluminal motion of PKS 0003-066 we measure 
corresponds to the middle to high end of the dis- 
tribution of superluminal motions in the sample of 
Kellermann et al. (2004, Fig.lO). 

3.1.2. The interaction between the moving and 
quasi-stationary components 

Also clear from the montage of images and Fig. 
9 is that the fast components C2 and C4 catch up 
to the slower component CI and interact strongly 
near epoch 1999. As this interaction occurs, the 
CI region clearly gets brighter. The component 
C3, closest to the core, appears stationary. 

Coexistence of faster and slower or quasi- 
stationary components are seen in many sources 
such as 4C 39.25 (Alberdi et al. 2000), Centau- 
rus A (Tingay et al. 2001), 3C 120 (Gomez et al. 
2001), 3C 279 (Jorstad et el. 2004) and 1803-^784 
(Britzen et al. 2005), which may be explained 
within relativistic time-dependent hydrodynamic 
models of "trailing shocks" (Agudo et al. 2001). 
In this scenario it is expected that the motion of 
inner components is slower than that of the outer 
primary component, as clearly seen in 3C 120 
(Gomez et al. 2001) and Centaurus A (Tingay et 
al. 2001). Our analysis of PKS 0003-066 suggests 
that C3 may be trailing shock for C2. 

Another explanation of the stationary fea- 
tures could be as bends of the jet trajectory in 
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a plane different to the observers plane, as sug- 
gested for 4C 39.25 by Alberdi et al. (2000). In 
PKS 0003—066 the component CI has a non-radial 
motion, from NE to SW, while brightning in total 
flux density. The component C3 is stationary in 
the distance from the core but not stationary in 
the direction perpendicular to the radial direction. 

The interaction between CI and C2 may be also 
evident in the variability of the polarization profile 
for CI. The VLBA polarization images at 15 GHz 
from 2003 February to 2005 June (Lister & Homan 
2005, MOJAVE database website) show that the 
component CI is significantly polarized, ~20% av- 
erage and up to 60% around the outer edge of the 
component, suggesting compression and enhance- 
ment of the local magnetic field at the interface 
with the ambient media. Such high fractional po- 
larization is not seen for the CI component in the 5 
GHz global VLBI image from 1995 May (Gabuzda, 
Pushkarev & Cawthorne 2000), with only an up- 
per limit of <11%, possibly because the compo- 
nent is not compressed enough at that stage. 

The brightning of a jet caused by internal 
shocks with a slightly varying jet speed is seen 
also in the well studied galactic X-ray binary jet 
SS 433. Mighari, Fender & Mcndcz (2002) show 
evidence for a hot continuum and Doppler-shifted 
iron emission lines from spatially resolved regions, 
suggesting that in situ reheating of the jet com- 
ponent takes place in a flow that moves with rel- 
ativistic bulk velocity more than 100 days after 
launch from the binary core of the SS 433 jet. 

3.1.3. Relation to large-scale structure 

As discussed earlier, the pc-scale VLBI jet kine- 
matics possibly supports jet bending. An interest- 
ing question is whether the structure on kpc-scales 
as seen with the VLA would be related to the 
pc-scale structure with an ^90° misalignment, as 
commonly seen in many compact core-dominant 
sources (e.g. Conway & Murphy 1993). In a 
report of a VLA calibrator position survey, Per- 
ley (1982) noted a single secondary component of 
PKS 0003-066 located at r = 1.7" (8.5 kpc), 
p.a=30°, but to our knowledge no kpc-scale image 
for this source has ever been published. We made 
two VLA images from archived data at 4.9 and 
8.5 GHz (Figs. 13(a) and (b)). We clearly detect 
the secondary component as pointed out by Perley 
(1982), although it appears very weak, approxi- 



mately 0.3 percent of the flux density from the 
unresolved core. The inner VLBI jet (< 0.8 mas) 
has an initial position angle centered at 7° with 
an oscillation of ±15°. On scales of 2-6 mas, the 
position angle is close to ~90°. On arcsec scales, 
the jet appears close to 0° once again. A natu- 
ral explanation of the extreme curvature of the jet 
is that we are looking at a part of a spiral-like 
trajectory, its appearance affected by projection 
onto the sky plane. An example of such orthogo- 
nal bending in a kpc-scale jet of a BL Lac object 
can be seen in 1803-1-784 as discussed by Tateyama 
et el. (2002), modeling a helical jet structure. 

4. Conclusion 

The main results of our multi-epoch VLBI im- 
age analysis of the BL Lac PKS 0003—066 are as 

follows: 

1. In contrast to previous studies of this source, 
VLBI components of the jet are highly variable, 
with proper motions of two components, C2 and 
C4, found to be as high as 0.6 mas/yr and 1.4 
mas/yr (~ 9 c and 22 c) respectively. The jet 
component CI is moving very slowly with ^0.1 
mas/yr or ~2.0 c. As the inner component C2 
approaches CI, the flux density of CI increases 3 
times, indicating an interaction between the com- 
ponents. 

2. A new disk-based recording system and 

software correlator for VLBI were verified by 
comparing imaging results to other established 
recorder /correlator systems. 
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U (10*^ A) 



Fig. 1. — {u,v) coverage of the global VLBI ob- 
servations of 2004 April 18 at 2.3 GHz. 



20 40 60 

UV radius (10^ A) 



Fig. 2. Visibility amplitude in Jy versus {u,v) 
distance, for the data shown in Figure 1. 
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-5 
Right Ascension (mas) 



Fig. 3. — Imago obtained from the global VLBI 
observations of 2004 April 18 at 2.3 GHz. The 
FWHM of the restored beam is 2.34 x 1.73 mas 
at position angle -13.9°, as indicated by the filled 
ellipse in the bottom left hand corner. Contours 
are drawn at -1.5, 1.5, 3, 6, 12, 24, 48, and 92% of 
1.31 Jy/beam, the peak flux density in the map. 






Fig. 4.— VLBA observations at 2.3 and 8.6 GHz 
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Fig. 8. — Time series mosaic of a selection of 8.6 
Fig. 7.— Time series mosaic of a selection of 2.3 GHz VLBI images of PKS 0003-066 (observing 

GHz VLBI images of PKS 0003-066. epochs are the same as 2.3 GHz). 




500 1000 1500 2000 2500 3000 3500 4000 
Modilied Julian Day +50,000 




500 1000 1500 2000 2500 3000 3500 4000 
Modified Julian Day +50,000 

Fig. 9. — (a) (top) Angular distance from the core 
of Gaussian components (CI, C2, C3 and C4) as 
a function of time (date from 50,000 Modified Ju- 
lian Date, i.e. 1995 October 10) for observations 
with the VLBA at 8.6 GHz and the global VLBI 
observations of 2004 April 18 at 2.3 GHz (MJD 
53,114). The lines are the linear least-square fits 
to the data, (b) (bottom) The same 8.6 GHz fit- 
ted lines as the above but, for comparison, with 
2.3 GHz components (triangles CI, circles C2, and 
square C4) and 15 GHz components shown by 
Kellermann et al. (2004) (filled squares compo- 
nent B, filled circles component C, and filled tri- 
angles component D). 




-2 -4 -6 -8 

Relative R.A. (mas) 

Fig. 10. — Position relative to the core of each 
component CI, C2, C3, C4 at all epochs derived 
from the observations with the VLBA at 8.6 GHz. 
Note that positional errors are estimated to be 
typically 0.2 and 0.5 mas for R.A. and Dec. re- 
spectively at 8.6 GHz taking one quarter of the 
beam sizes. 
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- CI 
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- C2 
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- C4 



•*ff,^;Tl'T-,-,-,-r,-,-4-,y, 
1000 1500 2000 2500 

Ivlodilied Julian Day +50.000 



.1 

3000 



3500 



Fig. 11. — Flux density variation for each VLBA 
8.6 GHz component as a function of time (date 
from 50,000 Modified Juhan Date, i.e. 1995 Octo- 
ber 10). The error bars on the flux density points 
are estimated at 10 percent of the flux density at 
each component. 
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500 1 000 1 500 2000 2500 3000 3500 4000 



Modified Julian Day +50,000 




500 1000 1500 2000 2500 3000 3500 4000 



Modified Julian Day +50,000 

Fig. 12.— (a) (top) The flux density of each VLBI 
component at 2.3 GHz. (b) (bottom) The total 
flux density of the VLBI core components and 
outer most components CI at 2.3 GHz and 8.6 
GHz. 




6000 4000 2000 -2000 -4000 -6000 

Right Ascension (mas) 
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4000 2000 -2000 -4000 

Right Ascension (mas) 
Map center: RA: 00 OS 13.893, Dec: -06 23 35.336 (2000.0) 



Fig. 13. — Kiloparsec scale structure of 
PKS 0003—066 imaged from VLA archive data: 

(a) (top) at 4.86 GHz on 1997 January 25, with 
the FWHM of the restored beam 1.04 x 5.26 (arc- 
sec) at P.A. -88.3°, and with contours at -0.075, 
0.075, 0.15, 0.3, 0.6, 1.2, 2.4, 4.8, 9.6, 19.2, 38.4, 
76.8% of 2.46 Jy/beam, the peak flux of the map, 

(b) (bottom) at 8.46 GHz on 2002 February 21, 
with the FWHM of the restored beam 2.99 x 1.94 
(arcsec) at P.A. -24.6°, and with contours at -0.05, 
0.05, 0.1, 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, 12.8, 25.6, 
51.2% of 2.76 Jy/beam, the peak flux of the map. 
Both observations were made with the VLA A- 
conflguration. 
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Table 1: Antennas participating in 2004 April 18 
global VLBI observations at 2.3 GHz 



Antenna 


Organisation*^ 


D (m) 


SEFD'' (Jy) 


Recording system 


ATCA 


ATNF 


5 X 22 


90 


LBADR 


Mopra 


ATNF 


22 


400 


LBADR 


Parkes 


ATNF 


64 


45 


LBADR 


Hobart 


U.Tas. 


26 


650 


Marks 


Ceduna 


U.Tas. 


30 


350 


LBADR 


Hartebeesthoek 


Hart.RAO 


26 


390 


Mark5 


Kashima 


NICT 


34 


400 


K5 



"ATNF is Australian Telescope National Facility. U.Tas. is University of Tasmania (Australia). Hart.RAO is Hartebeesthoek 
Radio Astronomy Observatory (South Africa). NICT is National Institute of Information and Communications Technology 
(Japan) 

''System Equivalent Flux Density 
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Table 2: Gaussian model of the image from global 
VLBI observations of 2004 April 18 at 2.3 GHz 



Component 




r(mas)'' 


PA(deg.)'' 


a (mas)*^ 






Core 


1.70 






1.24 


0.66 


74.1 


CI 


0.35 


6.70 


-82.8 


0.70 


0.00 


43.0 


C4 


0.32 


5.38 


-74.1 


1.68 


0.36 


49.3 



"Flux density in Janskys 

'r and PA arc the polar coordinates of the centre of the Gaussian relative to the core. Position angle is measured from north 
through east. 

'^a and b are the FWHM of the major and minor axes of the Gaussian. 
"^Position angle of the major axis measured from north through east. 
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Table 3: Parameters of the 2.3/8.6 GHz VLBI ob- 
servations from RRFID and VLBA archives 



Epoch 


Code 


Antenna 


2.3GHz beam 


8.6GHz beam 








(mas/mas/deg) 


(mas/mas/deg) 



1995 Oct 13 
1997 Jan 30 
1997 Mar 31 
1997 May 19 
1997 Jul 24 
1997 Sep 8 

1997 Dec 17 

1998 Feb 9 
1998 Apr 15 
1998 Jun 25 
1998 Aug 10 
1998 Dec 21 
2002 Jan 16 
2002 Jul 25 

2002 Sep 25 

2003 May 7 
2003 Jun 19 
2003 Sep 18 

2003 Dec 17 

2004 Mar 3 
2004 May 5 

2004 Jul 15 

2005 May 12 
2005 Jun 2 
2005 Jul 1 



RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RRFID 
RDV 34 
RDV 35 
RDV 38 
RDV 39 
RDV 41 
RDV 42 
RDV 43 
RDV 44 
RDV 45 
BP 118 
BP 118 
BP 118 



VLBA 

VLBA+6GA 

VLBA-I-6GA 

VLBA+6GA 

VLBA+6GA 

VLBA-I-6GA 

VLBA+7GA 

VLBA+7GA 

VLBA-I-7GA 

VLBA+7GA 

VLBA+7GA 

VLBA-I-7GA 

VLBA+8GA 

VLBA+8GA 

VLBA-I-8GA 

VLBA+9GA 

VLBA+8GA 

VLBA-I-8GA 

VLBA+9GA 

VLBA+8GA 

VLBA-I-8GA 

VLBA+8GA 

VLBA 

VLBA 

VLBA 



13.5/3. 97/-15.3 

8.02/2.88/-4.76 

6.55/2.85/-0.75 

8.89/3. 20/-6. 20 

8.41/3.26/-3.10 

11.0/3.25/-11.5 

9.59/2. 83/-13.7 

10.2/2.98/-13.5 

7.45/3.01/-2.55 

10.2/3.28/-3.67 

8.98/2.49/-8.97 

8.60/3.57/1.04 

9.20/2.46/-5.09 

6.80/4.45/26.4 

6.82/1.91/1.73 

2.64/1.83/32.0 

7.89/4.68/26.5 

5.61/2.24/-0.761 

5.51/2.41/6.17 

5.96/2.09/-1.33 

6.96/2.80/-11.8 

5.50/2.08/-4.83 

6.33/2.36/5.05 

6.45/2.49/-0.392 

6.24/2.55/-0.351 



2.29/0.955/-1.13 

2.03/0.745/-5.85 

1.71/0.783/1.67 

2.31/0.849/-5.68 

2.24/0.837/-3.52 

2.29/0.867/-11.6 

2.18/0.718/-10.7 

2.47/0.742/-11.4 

2.04/0.799/-4.24 

2.59/0.712/-3.90 

1.73/0.653/-4.29 

2.35/0.777/-0.392 

2.40/0. 598/-6. 26 

2.44/1.50/-11.0 

1.85/0.500/1.45 

0.841/0.554/11.3 

2.01/1.27/20.8 

1.67/0.616/-1.48 

0.808/0.530/51.9 

1.69/0.580/-1.38 

1.92/0.787/-12.3 

1.54/0.564/-7.10 

1.75/0.636/2.99 

1.82/0.658/-2.16 

1.65/0.720/0.581 



Note. — VLBA consists of 10 antennas of 25 m. GA stands for geodetic antennas as described in Tateyama & Kingham (2004); 
6GA = Gc, Gn, Kk, Mc, On, and Wf: 7GA=6GA+Ny; 8GA=7GA+Ts-Wz-Gn: 9GA= 8GA+Tc. Gc is Gilcreek (Alaska, 26m), 
Gn is NRAO20 (Grecnbank, 20m), Kk is Kokee (Hawaii, 20m), Mc is Medicina (Italy, 32 m), Ny is Ny Alesund (Norway, 20 
m). On is Onsala (Sweden, 20 rn), Tc is Tigo (Chile, 6m), Ts is Tsukuba (Japan, 32 m), Wf is Westford (Massachusetts, 18 
m), and Wz is Wettzell (Germany, 20 m). 
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Table 4 

Gaussian models for ial\ges of the VLBI observations from RRFID and VLBA archives 



Epoch (MJD) 


GHz 


Gomp.ID 


5(Jy) 


r(nias) 


PA(deg.) 


a (mas) 


b/a 


tpideg.) 


lyyo uct io (^ouuuo J 


Z.O 


Core 


i.oOi 






\J.t<0 


(~i nn 
U.UU 


7^? 9 
/O.Z 






L^z 


U.DOo 


z.UU 


-Di.4 


i.oi 


n n9 

u.yz 


O.y 






v^i 


u.zyu 


0.4/ 


- ( O.O 


1 f;9 
i.OZ 


n QQ 

U.OO 


/I 7 /I 
4 ( .4 




o.O 


Core 


1 A 

i.044 






U.OO 


n /1 9 
U.4z 


9/^ Px 
-ZU.D 






Oo 


u.oy 


i.U / 


-00.4 


1 70 

i. (^z 


n QQ 
U.OO 


O 

00. Z 






CA 




O.Z / 


- ( Z.O 


Q f: /I 
o.o4 


n Q 1 
U.oi 


Q 

00. 


lijiji Jan oU yoiJ^iv) 


Z.O 


Gore 


i.4oo 






1 9^^ 
i.ZD 


n 79 


(^9 9 
-Dz.z 






l^Z 


u.oyo 


z.zo 


- / 4.D 


1 09 

i.yz 


fi (1(^ 
U.UU 


111 








U.Z ( 


o.iy 


7C^ 1 
- ( O.i 


9 C^/1 
Z.04 


n 71 
U. /i 


9/1 1 
z4. i 




O.t) 


Gore 








U.04 


n 1 Q 
U.io 


1 O (? 

ly.o 






L^O 


u.iyy 


U.oO 


A^ P. 
-4/.D 


n /( c 
U.4o 


n 7Q 

U. ( 


/I Q /) 
40.4 








U.zoz 


9 11 
Z.ii 


7ri 
- ( U.O 


1 cn 


n 

U.OU 


7Q A 
- 1 0.4 








U.UoO 


0. (>) 


7Q n 
- / o.y 


O /I Q 

z.4o 


n QQ 
U.OO 


1 1 n 
ii.U 


lyy/ Mar oi (oUooyj 


Z.O 


Gore 


1.661 






1 07 
i.Z (^ 


n ^^^^ 
U.OO 


-Oz.i 






*^z 


U.Oi ( 


9Q 
z.Zo 


7/1 Q 
- ( 4.0 


1 OQ 

i.yo 


n 7Q 

u. / o 


1 n /i 

-iU.4 








U.zz / 


^ c;i 


7c: Q 
- (O.O 


9 17 
Z.i ( 


n t^7 
U.O ( 


9c; n 
zO.U 




o.D 


Gore 


1 OQQ 






U.OO 


n QQ 

U.OO 


Q Q 
0.0 








n 01 c:^ 


U.oy 


/1A Q 
-40. o 


U.04 


n c:q 
U.Oo 


-00.0 






L^z 


n OOQ 


9 Q9 
Z.OZ 


7n 9 

- (U.Z 


1 

i.Oo 


n /I Q 
U.4o 


-DD.Z 








U.U ( o 


D.UU 


7Q O 

- ( O.y 


9 Q7 
Z.O ( 


n ^^n 

U.oy 


97 Q 


lyy/ iviay ly (^ouoooj 


O Q 
Z.O 


Gore 


i.oo / 






1 /i A 
1.40 


n t^n 
U.OU 


^O Q 

-Oz.o 






v^z 




9 1^1 


7 A Q 
- ( 4.y 


1 .UU 


n nn 

U.UU 


-0.0 








U.zoo 


o.zo 


7^ O 


Q 1 (^ 
O.iO 


n /1 9 
U.4z 


/I (^ 
4.0 




O.D 


Gore 


1 1/10 






U.oy 


n Q /( 
U.o4 


1 f; n 

io.y 






Oo 


n oi 7 
U.zi ( 


U.OO 


/in n 


U.OO 


n 7/1 


fiA A 
-04.4 






L^z 


U.zzo 


9 Q /I 


7n Q 

- (U.O 


9 f~19 

z.Uz 


n f^i 
U.Oi 


7/1 Q 
- (^ 4.0 






Ui 


U.Uol 


o.y / 


7*5 Q 
- ( O.O 


9 1 n 
z.iU 


n Qn 
U.oU 


Q7 


1997 Jul 24 (50654) 


2.3 


Core 






0.48 


1.45 


0.42 


-60.0 






G2 


0.375 


2.58 


-73.4 


1.63 


0.00 


-1.5 






Gl 


0.252 


5.34 


-75.3 


3.10 


0.35 


10.4 




8.6 


Gore 


1.323 






0.56 


0.37 


20.4 






G3 


0.190 


0.90 


-51.0 


0.46 


0.07 


-77.6 






G2 


0.234 


2.40 


-71.6 


1.99 


0.62 


-79.9 






Gl 


0.081 


5.93 


-74.7 


2.19 


0.74 


23.9 


1997 Sep 8 (50700) 


2.3 


Gore 


1.56 






1.84 


0.69 


-28.8 






G2 


0.408 


2.74 


-73.9 


2.15 


0.06 


13.7 






Gl 


0.221 


5.59 


-74.1 


3.28 


0.00 


14.3 




8.6 


Gore 


1.426 






0.62 


0.36 


19.6 






G3 


0.201 


0.91 


-50.2 


0.49 


0.00 


84.9 






G2 


0.252 


2.44 


-70.8 


2.23 


0.61 


-74.7 






Gl 


0.081 


6.09 


-73.8 


2.47 


0.53 


16.7 
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Table 4 -Continued 



Epoch (MJD) GHz Comp.ID S'(Jy) r(mas) PA(deg.) a(mas) h/a ifjideg.) 



1997 Dec 17 (50800) 2.3 Core 1.538 

C2 0.404 

CI 0.251 

8.6 Core 1.390 

C3 0.375 

C2 0.209 

CI 0.085 

1998 Feb 9 (50854) 2.3 Core 1.358 

C2 0.451 

CI 0.189 

8.6 Core 1.559 

C3 0.282 

C2 0.200 

CI 0.071 

1998 Apr 15 (50919) 2.3 Core 1.386 

C2 0.428 

CI 0.258 

8.6 Core 1.731 

C3 0.272 

C2 0.208 

CI 0.070 

1998 Jun 25 (50989) 2.3 Core 1.526 

C2 0.347 

CI 0.176 

8.6 Core 1.714 

C3 0.427 

C2 0.167 

CI 0.057 

CO? 0.016 

1998 Aug 10 (51036) 2.3 Core 1.463 

C2 0.359 

CI 0.250 

8.6 Core 1.863 

C3 0.296 

C2 0.208 

CI 0.074 

1998 Dec 21 (51169) 2.3 Core 1.437 

C2 0.344 

CI 0.202 

8.6 Core 2.160 

C3 0.428 







1.56 


0.46 


-58.3 


2.65 


-75.8 


2.14 


0.13 


1.0 


5.38 


-75.6 


1.94 


0.75 


17.0 






0.65 


0.00 


22.6 


0.60 


-53.2 


1.05 


0.41 


-59.7 


2.63 


-73.7 


2.03 


0.63 


-79.2 


6.07 


-76.6 


2.17 


0.72 


53.7 






1.46 


0.65 


-49.2 


2.75 


-74.9 


2.08 


0.62 


-0.2 


5.70 


-75.6 


2.16 


0.00 


18.7 






0.63 


0.08 


22.7 


0.73 


-60.6 


0.90 


0.43 


-50.8 


2.75 


-72.7 


2.14 


0.69 


-80.6 


6.22 


-76.5 


1.59 


0.94 


12.4 


_ 


_ 


1.19 


0.66 


-67.9 


2.57 


-74.4 


1.86 


0.46 


22.1 


5.39 


-74.4 


2.18 


0.63 


53.8 






0.62 


0.13 


20.6 


0.79 


-64.9 


0.88 


0.53 


-30.4 


2.77 


-74.3 


2.26 


0.63 


-80.6 


6.24 


-75.2 


1.90 


0.58 


22.5 






1.84 


0.79 


-47.7 


3.30 


-68.3 


2.74 


0.38 


1.5 


5.73 


-79.4 


4.93 


0.00 


-3.8 






0.61 


0.00 


21.4 


0.59 


-72.0 


1.33 


0.25 


-59.8 


2.86 


-75.2 


2.11 


0.56 


-79.6 


6.21 


-76.7 


1.95 


0.47 


32.2 


9.67 


-76.4 


1.00 


0.00 


26.4 






1.27 


0.41 


-79.7 


2.69 


-74.2 


1.55 


0.00 


28.2 


5.34 


-75.4 


2.11 


0.51 


56.7 






0.59 


0.20 


19.5 


0.86 


-69.4 


0.93 


0.53 


-25.8 


2.84 


-75.3 


2.40 


0.53 


-88.3 


6.30 


-77.0 


1.47 


0.95 


64.1 






1.69 


0.50 


86.2 


2.95 


-71.6 


1.51 


0.54 


52.0 


5.70 


-76.5 


2.83 


0.29 


29.0 






0.56 


0.51 


15.8 


0.86 


-68.5 


1.68 


0.43 


-78.9 
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Table 4— Continued 



Epoch (MJD) GHz Comp.ID S{3y) r(mas) PA(deg.) o(mas) b/a tl>{deg.) 



C2 0.086 

CI 0.088 

2002 Jan 16 (52291) 2.3 Core 1.211 

C2 0.193 

CI 0.503 

8.6 Core 1.202 

C4 0.197 

C3? 0.015 

C2 0.036 

CI 0.203 

2002 Jul 25 (52481) 2.3 Core 1.000 

CI 0.463 

8.6 Core 1.979 

C4 0.133 

C2 0.072 

CI 0.220 

2002 Sep 25 (52543) 2.3 Core 1.201 

C4 0.390 

CI 0.507 

8.6 Core 1.794 

C3 0.286 

C2 0.079 

CI 0.218 

2003 May 7 (52767) 2.3 Core 1.384 

C4 0.566 

CI 0.766 

8.6 Core 2.088 

C4 0.176 

C3 0.122 

CI 0.291 

2003 Jun 19 (52809) 2.3 Core 2.362 

C4 0.688 

8.6 Core 2.113 

C4 0.181 

CI 0.401 

2003 Sep 18 (52886) 2.3 Core 1.108 

C4 0.602 

CI 0.702 

8.6 Core 2.027 

C2 0.209 

C3 0.408 



3.42 


-80.3 


1.60 


0.00 


17.8 


6.30 


-77.7 


4.17 


0.04 


24.0 






1.72 


0.86 


-43.7 


3.19 


-66.0 


1.70 


0.00 


4.4 


5.87 


-80.4 


2.09 


0.65 


-3.5 






0.80 


0.32 


5.5 


1.48 


-59.8 


1.33 


0.73 


5.8 


3.32 


-70.5 


3.15 


0.00 


4.1 


4.26 


-83.4 


2.94 


0.00 


76.4 


6.53 


-77.0 


2.06 


0.42 


1.9 






1.82 


0.56 


-76.4 


5.60 


-79.7 


2.38 


0.69 


82.2 






0.86 


0.56 


13.6 


2.35 


-66.9 


1.35 


0.47 


16.3 


5.04 


-70.0 


1.79 


0.00 


21.8 


6.61 


-78.3 


1.58 


0.36 


12.1 






1.45 


0.23 


-15.2 


2.43 


-79.2 


3.73 


0.25 


-50.8 


6.14 


-80.9 


1.79 


0.00 


-41.1 






0.81 


0.30 


-0.1 


1.29 


-80.9 


2.51 


0.56 


-6.2 


5.01 


-78.3 


3.54 


0.38 


17.8 


6.69 


-79.5 


1.49 


0.53 


13.0 






1.97 


0.43 


1.5 


1.45 


-74.6 


4.63 


0.57 


-27.2 


6.14 


-81.2 


2.48 


0.56 


7.6 






0.83 


0.36 


-0.1 


2.40 


-69.5 


4.14 


0.00 


85.9 


1.04 


-94.7 


0.46 


0.00 


89.2 


6.70 


-78.8 


1.41 


0.74 


14.3 






3.77 


0.57 


-72.4 


6.56 


-78.7 


7.95 


0.00 


8.1 






0.87 


0.47 


-8.2 


3.07 


-79.0 


3.37 


0.00 


-70.5 


6.47 


-81.3 


1.79 


0.69 


-4.6 






2.06 


0.25 


-20.6 


1.66 


-80.5 


5.13 


0.52 


-43.0 


6.34 


-81.0 


2.79 


0.46 


-17.5 






0.74 


0.43 


-1.8 


6.22 


-73.9 


3.19 


0.00 


-72.5 


1.797 


-94.8 


3.26 


0.00 


-50.9 
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Table 4 - Continued 



Epoch (MJD) GHz Comp.ID S'(Jy) r(mas) PA(deg.) a(mas) b/a i^ideg.) 



CI 0.186 

2003 Dec 17 (52991) 2.3 Core 1.199 

C4 0.295 

CI 0.612 

8.6 Core 2.229 

C4 0.060 

C3 0.185 

C2 0.058 

CI 0.249 

2004 Mar 3 (53067) 2.3 Core 1.068 

C4 0.286 

CI 0.559 

8.6 Core 1.768 

C3 0.147 

CI 0.251 

2004 May 5 (53131) 2.3 Core 1.829 

C4 0.128 

CI 0.883 

8.6 Core 2.364 

C4 0.171 

C3 0.109 

CI 0.349 

2004 Jul 15 (53201) 2.3 Core 1.351 

C4 0.352 

CI 0.694 

8.6 Core 2.481 

C4 0.128 

C5? 0.356 

C3 0.188 

CI 0.332 

2005 May 12 (53493) 2.3 Core 1.473 

C4 0.121 

CI 0.641 

8.6 Core 2.098 

C4 0.201 

C3 0.467 

CI 0.167 

2005 Jun 2 (53523) 2.3 Core 1.602 

C4 0.354 

CI 0.748 

8.6 Core 2.268 



6.78 


-82.2 


1.41 


0.49 


-30.7 


_ 


_ 


1.18 


0.00 


-47.9 


2.61 


-70.7 


3.01 


0.64 


-34.7 


6.38 


-80.8 


2.04 


0.53 


-17.0 






0.62 


0.53 


8.7 


3.03 


-74.2 


2.39 


0.27 


40.3 


0.76 


-61.9 


1.81 


0.00 


40.7 


6.63 


-70.2 


0.59 


0.00 


89.5 


6.77 


-82.0 


1.42 


0.54 


-26.3 






1.37 


0.00 


-27.1 


2.24 


-69.9 


3.92 


0.12 


-60.0 


6.41 


-81.0 


1.64 


0.55 


-30.0 


_ 


_ 


0.49 


0.46 


26.8 


0.71 


-40.2 


0.75 


0.00 


46.5 


6.87 


-79.4 


1.59 


0.59 


6.6 






1.86 


0.55 


-25.4 


3.10 


-71.3 


13.30 


0.00 


-18.6 


6.07 


-82.9 


2.70 


0.66 


-12.9 






0.51 


0.00 


23.1 


4.53 


-78.8 


6.80 


0.42 


-37.5 


1.56 


-84.8 


2.01 


0.22 


-19.0 


6.80 


-82.4 


1.45 


0.59 


13.2 


_ 


_ 


1.14 


0.70 


-9.7 


2.29 


-73.7 


4.00 


0.52 


-49.0 


6.48 


-81.2 


2.08 


0.61 


2.9 






0.52 


0.44 


21.0 


4.74 


-78.7 


3.68 


0.92 


-6.7 


0.59 


-56.4 


0.61 


0.00 


53.4 


1.16 


-97.2 


3.15 


0.11 


-39.0 


6.89 


-80.7 


1.43 


0.70 


1.4 






1.28 


0.81 


-61.9 


3.01 


-70.1 


0.42 


1.00 


-14.4 


6.29 


-81.0 


2.58 


0.46 


3.7 






0.76 


0.32 


9.65 


5.69 


-85.0 


6.48 


0.49 


63.5 


0.91 


-103.3 


2.15 


0.42 


27.3 


7.02 


-81.4 


1.58 


0.00 


-51.0 






1.23 


0.85 


9.5 


2.42 


-71.7 


3.06 


1.00 


-90.0 


6.44 


-80.9 


2.37 


0.49 


9.0 






0.81 


0.38 


2.9 
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Table 4 — Continued 



Epoch (MJD) 


GHz 


Comp.ID 




r(mas) 


PA(deg.) 


a(mas) 


b/a 


■0(deg.) 






C4 


0.069 


5.87 


-92.1 


4.05 


0.00 


34.6 






C3 


0.269 


0.69 


-65.1 


0.56 


0.34 


87.6 






CI 


0.277 


6.88 


-80.7 


1.62 


0.59 


-2.5 


2005 Jul 1 (53554) 


2.3 


Core 


1.752 






1.27 


0.84 


-87.8 






C4 


0.216 


3.01 


-68.4 


1.48 


1.00 


-99.6 






CI 


0.765 


6.34 


-81.1 


2.08 


0.64 


16.0 




8.6 


Core 


2.262 






0.79 


0.42 


4.1 






C3 


0.181 


0.73 


-83.2 


1.90 


0.00 


3.2 






CI 


0.367 


6.78 


-81.0 


1.98 


0.57 


6.3 



23 



Tabic 5: Parameters of the 2.3/8.6 GHz band 
VLBI maps shown in Fig. 3-5. 



2.3GHz peak 2.3GHz contour 8.6GHz peak 8.6GHz contour 



Epoch 


(Jy/beam) 


% of peak 


(Jy/beam) 


% of peak 


2002 Jul 25 


1.07 


-2 2 4 8 16 32 64 


1.77 


-0.4 0.4 0.8 1.6 3.2 6.4 12.8 25.6 


2002 Sep 25 


1.18 


-1 1 2 4 8 16 32 64 


1.53 


-1 1 2 4 8 16 32 64 


2003 May 7 


1.09 


-1.5 1.5 3 6 12 24 48 96 


1.26 


-1 1 2 4 8 16 32 64 


2003 Jun 19 


1.8 


-1 1 2 4 8 16 32 64 


1.93 


-1.5 1.5 3 6 12 24 48 96 


2003 Sep 18 


1.16 


-1.5 1.5 3 6 12 24 48 96 


1.64 


-0.75 0.75 1.5 3 6 12 24 48 96 


2003 Dec 17 


1.22 


-1 1 2 4 8 16 32 64 


1.03 


-1.5 1.5 3 6 12 24 48 96 


2004 Mar 3 


1.09 


-0.75 0.75 1.5 3 6 12 24 48 96 


1.55 


-0.5 0.5 1 2 4 8 16 32 64 


2004 May 5 


1.66 


-0.5 0.5 1 2 4 8 16 32 64 


2.51 


-0.5 0.5 1 2 4 8 16 32 64 


2004 Jul 15 


1.29 


-0.75 0.75 1.5 3 6 1 12 24 48 96 


1.87 


-0.5 0.5 1 2 4 8 16 32 64 


2005 May 12 


1.56 


-1 1 2 4 8 16 32 64 


1.83 


-0.75 0.75 1.5 3 6 12 24 48 96 


2005 Jun 2 


1.56 


-1 1 2 4 8 16 32 64 


1.96 


-0.5 0.5 1 2 4 8 16 32 64 


2005 Jul 1 


1.59 


-0.5 0.5 1 2 4 8 16 32 64 


1.90 


-0.5 0.5 1 2 4 8 16 32 64 
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Table 6: Proper motions modeled from multi 
epoch VLBI data at 8.6 GHz 



Component 


to(Modified J-D.)"* 


to(Epoch)'^ 


/u(mas/yr) 


l^app 


CI 


32470 ± 1866 


1947.79 ± 5.11 


0.12 ± 0.01 


2.0 ± 0.2 


C2 


49185 ± 75.3 


1993.54 ± 0.21 


0.58 ± 0.02 


9.3 ± 0.3 


C3 


24942 ± 26268 


1927.16 ± 71.97 


0.01 ± 0.01 


0.2 ± 0.2 


C4 


51977 ± 42.9 


2001.18 ± 0.12 


1.35 ± 0.05 


21.6 ± 0.8 



"Zero Separation Time 
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